Weak First-Order Superfluid— Solid Quantum Phase Transitions 
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We study superfluid-solid zero-temperature transitions in two-dimensional lattice boson/spin 
models by Worm-Algorithm Monte Carlo simulations. We observe that such transitions are typ- 
ically first-order with the exception of special high-symmetry points which require fine tuning in 
the Hamiltonian parameter space. We present evidence that the superfluid-checkerboard solid and 
superfluid-valence-bond solid transitions at half-integer filling factor are extremely weak first-order 
transitions and in small systems they may be confused with continuous or high-symmetry points. 
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Recently, there has been an increased interest in 
superfluid-solid (SF-S) quantum phase transitions in lat- 
tice boson/spin systems [l|, H; IM Q- [By 'solid' we 
mean an insulating state featuring the broken transla- 
tion symmetry — like checkerboard solid/ antiferromagnet 
(CB) or valence bond solid (VBS), as opposed to the 
Mott insulator which preserves the translation symme- 
try.] On one hand, this interest is stimulated by experi- 
mental pers^pectives of studying such transitions in opti- 
cal lattices Q, on the other hand, the SF-VBS transition 
in a (2 -I- l)-dimensional system has been argued to be 
the simplest example of qualitatively new type of quan- 
tum criticality, that does not fit the standard Landau- 
Ginzburg- Wilson paradigm [j|. Intriguing data on the 
SF-S transitions were obtained by direct Monte Carlo 
simulations of quantum systems 1, 2]. It was observed l] 
that the SF-CB transition in the hard-core bosonic model 
with nearest- and next-nearest-neighbor interaction re- 
mains remarkably insensitive to the explicit Heisenberg 
symmetry breaking. The transition is numerically indis- 
tinguishable from that of the Heisenberg spin- 1/2 model: 
a degenerate (hysteresisless) transition. In simulations of 
the 2D quantum XY model with ring exchange 2], the 
SF-VBS transition was interpreted as the second-order 
one, which suggested its understanding in terms of the 
deconfined quantum critical point |j| . 

In this Letter, we perform a careful finite-size analy- 
sis of the SF-CB and SF-VBS transitions by simulating a 
generalized (2-|-l)-dimensional J-current model 5], which 
is a discrete-imaginary-time analog of the quantum bo- 
son/spin lattice system. The simplicity and flexibility of 
our model in combination with the efficient Worm Algo- 
rithm allow us to arrive at a definitive conclusion that 
in SF-CB and SF-VBS cases we are dealing with anoma- 
lously weak first-order phase transitions. Moreover, the 
two transitions are remarkably similar to each other. In 
both cases, small enough critical systems mimic the be- 
havior of a highly symmetric model with broken symme- 
try. The SF-CB case corresponds to the 0(3)-symmetry 
of the Heisenberg model, while in the SF-VBS case we 
clearly see a quasi-0(4) behavior that manifests itself as 
a coexistence at the critical point of the superfluid re- 
sponse and both (x- and y-) VBS orders, in arbitrary 



proportions. 

We confine ourselves to the case when the statistics 
of (2-1-1 )-dimensional particle trajectories in imaginary 
time (worldline configurations) is positive definite. What 
groundstates can emerge under this condition? The state 
with chaotic (unstructured) typical worldline configura- 
tion is SF. Indeed, the absence of structure implies fiuc- 
tuations of the worldline winding numbers, W , and thus 
a finite superfluid response which is given by the mean 
square of VF Q, see Eq. ^ below. [One can hardly 
extend this argument to cases when the configuration 
weight is not positive definite, since the notion of a typ- 
ical configuration becomes vague.] In SF, the U(l) sym- 
metry is broken (at least in a topological sense). The only 
way to restore this symmetry is to suppress fluctuations 
of W which seems to be impossible without structur- 
ing worldline configurations in such a way that for each 
imaginary-time moment the position of each worldline 
in the corresponding real-space plane can be unambigu- 
ously associated at the microscopic level with one of the 
lattice sites/bonds, and vice versa. If the total num- 
ber of the worldlines is not equal to a multiple of the 
number of sites/bonds, the structured worldline configu- 
ration immediately implies a broken translation symme- 
try. This consideration leads to the conjecture that for 
models with the positive definite statistics of worldline 
configurations and non-integer filling factor the generic 
groundstate should feature an order, either SF or solid, 
or both (supersolid) . Groundstates with none symmetry 
being broken ( "quantum disorder" ) may then occur only 
as critical points separating the ordered states. 

In the path-integral representation, we see no quali- 
tative difference between the site- and bond-based solids 
since both are characterized by the worldline structuring, 
in the above-mentioned sense. In either case, zero-point 
fiuctuations necessarily include permutations of two and 
more lines, and thus on large scales such microscopic de- 
tails as the most probable positioning — on sites or on 
bonds — of the worldlines can hardly be relevant to the 
universal properties of the SF-S transition. The only 
property that seem to be crucial is the symmetry of the 
emerging lattice. 

In terms of algorithmic simplicity and numerical ef- 
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FIG. 1: SF stiffness ps and VBS order parameter B depen- 
dence on the coupling strength D for L = 16 (squares), L = 32 
(open circles on inset) and L = 64 (filled circles on inset) sys- 
tems. A closer look at the transition point is provided in the 
insets. Error bars are shown for all points (in some cases these 
are smaller than the symbol size). 



ficiency, classical id -f l)-dimensional analogs of d- 
dimensional quantum systems offer a significant advan- 
tage. This approach was successfully used previously in 
the studies of disordered 0, Qi 111 and two-component 
systems !9| . In addition, classical models offer more free- 
dom in "designing" effective models with complex phase 
diagrams. The so-called J-current model of Ref. ^| is 
obtained by considering trajectories in discrete imagi- 
nary time. Let n = (x, r) denote points of the (d -f 
l)-dimensional space-time lattice, and integer currents 
J^(n) with V = x,y^f specify how many particles are 
going from site n in direction v. In this language, cur- 
rents in the time-direction represent occupation num- 
bers, and currents in the space directions represent hop- 
ping events. The continuity of trajectories requires that 
E.[J.(n)-J.(n-j.)] = 0. 

The simplest J-current model at half-integer filling fac- 
tor is obtained by writing the potential energy term in 
the particle-hole symmetric form 



Sj^2J J2 (jrin) - I) (jfin + i^)-l 



n.y^T 



(1) 



and restricting currents along f-bonds to take on just two 
values, and 1. The kinetic energy term is simply 



Sk^K Y, \Un)\, 



(2) 



with the restriction that allowed values for spatial cur- 
rents are 1, 0, and —1 values. To exclude somewhat 
pathological cases with two hopping events happening 
at the same space-time point, we further require that 
X^i^T^f l'^i'('^)l + \Jv{t^ — i^)\ < 1- Finally, we introduce 
interactions between the spatial currents on n.n. bonds 



which favor VBS 
Sd = -DY, l^.(n)|^(|J.(n + /i)| + |J.(n-/i)|). 

(3) 
Equal-time coupling [n ^ t) favors simultaneous hopping 
events of particles on the same plaquette and is reminis- 
cent of the ring exchange term in quantum models j^]. 
Phonon mediated exchange is another known mechanism 
of dimerization in spin Pierles systems [Iflj , and in Eq. 
it is represented by coupling between bonds connecting 
the same sites and shifted in the time direction. With all 
three terms combined, S = Sj + Sk + Sd , the resulting 
model has SF, CB, and VBS states in its phase diagram. 
First we study the SF-VBS transition along the J = 0, 
K = 0.4 line. Superfluid stiffness is determined by the 
statistics of winding number fluctuations |6j 



{W^)/2L 



(4) 



and the VBS order parameter is characterized by the 
staggered distribution of spatial currents along x and y 
directions 



B, = L-f'^+i) 






\Mn)\e^ 



(5) 



where v — x and t^ = ^ for q = (tt, 0, 0) and q = (0, n, 0), 
respectively. We find it convenient to introduce a sin- 
gle VBS order parameter with positive definite estimator 
which takes on finite value ~ 0{1) in the VBS phase, 
B = \Bx\ + \By\. For completeness, we also define here 



the CB order parameter as 

M,(q=(^,7r,0)) = L-('^^ 



'Jf{n)e' 



(6) 



and M = \Mr\. 

In Fig. n] we show rescaled data for the SF stiffness 
Ps/ps{D = 0) and VBS order parameter. The main plot 
for L = 16 demonstrates strong suppression of ps and 
B near the critical point [Dc « 0.5705(2)] which is typ- 
ical for continuous phase transitions. Similar behavior 
was reported previously for the ring-exchange model in 
Ref. |3]. However, in the insets we clearly see that finite- 
size scaling is incompatible with the second-order transi- 
tion scenario — the curves Ps{D) for different sizes L inter- 
sect each other without any further rescaling indicating 
that large systems are more ordered in the vicinity of the 
critical point. Similar behavior is observed for the curves 
B{D). The most obvious scenario is then a first-order 
transition where the intersection of finite-size curves at 
the critical point is allowed. Apparently, the transition 
is weakly first-order because (i) both order parameters 
are strongly suppressed at Dc^ and (ii) simulations for 
V^ = 32"^ system do not show any hysteresis, though the 
autocorrelation time is very long at D^- 

The other surprising fact is that the region where ps for 
L = 16 is below the corresponding curves for L = 32 and 
L = 64 is rather extended, while in first-order transitions 
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FIG. 2: Statistics of fluctuating SF and VBS order parameters. Left panel: points of the "black square" are obtained for MC 
configurations separated by equal long time intervals for L — 16 and D = 0.5727. Middle panel: SF stiffness is calculated as a 
function of B for various system sizes. From top to bottom curves are plotted at corresponding finite-L "transition points" (see 
the text), for (L = 8, L> = 0.5742), (L = 12, D = 0.5738), (L = 16, L» = 0.5727), (L = 24, D = 0.5715), (L = 32, D = 0.5712), 



(L = 48, D — 0.5705). Right panel: the coarse-grained distribution of the average density of points in the B^ 

the B = \Bi\ + |By|=const lines for the same set of system sizes and values of D as in the middle panel. 
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it is expected to be macroscopically small. It appears as 
if the superfluid order parameter ^ experiences anoma- 
lously large fluctuations in small systems. To explain 
it we first speculate (and later prove numerically) that 
in the vicinity of the critical point the system is best 
described by the four-dimensional order parameter, S, 
and the 0(4)-symmetry is broken at Dc ■ Formally, B^, 
By and two components of ^ represent observable (lin- 
early independent) projections of S; correspondingly, in 
the {Bx, By, '^)-space the 0(4)-sphere is seen as a four- 
dimensional surface with the sphere topology. In this 
scenario, if the 0(4)-symmetry is exact then any value of 
S is equally probable at D^ i.e. solid orders along both 
spatial direction and superfluidity may coexist. 

An analogous well-known example of the 0(3)- 
symmetric point is provided by the SF-CB transition 
in the 2D quantum XXZ-antifcrromagnet with n.n. ex- 
change interactions. In this case, the XY order parameter 
Sx + iSy = \1/ and the CB order parameter Sz = M form 
a three-dimensional vector S. The critical point itself is 
described by the SU(2)-symmetric Heisenberg Hamilto- 
nian. Since the 0(3)-symmetry is spontaneously broken 
in the ground state of the Heisenberg model, at the tran- 
sition both ps and M change discontinuously in the ther- 
modynamic limit, but this change occurs without energy 
barriers and is preceded by anomalously large fluctua- 
tions and long autocorrelation times in finite systems. 

If the outlined picture is correct, then small pertur- 
bations which explicitly break the 0(n) symmetry at the 
microscopic level should result in generic weak first-order 
transitions. Indeed, in the spontaneously ordered state 
all renormalizations are finite. Thus symmetry breaking 
perturbations, which couple to the long-range order, re- 
sult in the non-flat macroscopic energy/effective action 
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FIG. 3: System size dependence of ps{B — Q)/pa[D = 0) 
(circles) and Bmax (squares) from Fig. 



profile for the order parameter. Phase transitions are, 
then, between energy minima separated by macroscopic 
barriers which may, however, remain relatively weak in 
small systems. The J-current model studied here has no 
microscopic symmetries mandating exact 0(4)- or 0(3)- 
symmetry of the critical point. We conclude then that 
SF-VBS and SF-CB transitions are expected to be first 
order, and in the rest of the paper we present evidence 
that this is indeed so. 

First, we demonstrate that relatively "small" (hun- 
dreds of particles!) systems in the vicinity of the critical 
point behave as if {B^, By, Vf) fluctuations are confined to 
some surface, not volume, and it is not possible to have 
all three order parameters fluctuating to zero. In the 
left panel of Fig.|21we show statistical fluctuations of the 
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FIG. 4: The coarse-grained distribution of the CB order pa- 
rameter for various system sizes; from right to left: (L = 
8, J = 0.450), (L = 16, J = 0.4205), (L = 24, J = 0.4189), 
(L = 32, J = 0.4186). 



VBS order parameters for L — 16 and D = 0.5727. The 
density of points in the Bx,By "black square" is nearly 
homogeneous, i.e. the system is equally likely to be found 
with small or relatively large VBS order oriented at any 
angle relative to x, y-directions. Moreover, the boundary 
shape suggests that Bx and By fluctuations happen on 
the B — /(l^l) surface. The middle panel of Fig. |21 shows 
that indeed the values of B and ps are strongly correlated 
[ps{B) is defined as the average, see Eq. (0J, over config- 
urations with the VBS order parameter e {B — 5,B + 5)]. 
We see that, Ps is largest when B ^ 0, i.e. despite large 
fluctuations of all order parameters, they never vanish 
simultaneously. 

Finally, in the right panel we plot the average den- 
sity of points in the i?:^, By-plane along the i?=const 
cuts for various system sizes. In each case, wc adjusted 
the coupling constant D so that the distribution function 
P{B) is "maximally flat" , or the two maxima are at equal 
heights. Normalization was set to have the large-B max- 
imum equal unity. The density proflles are nearly flat 
for L = 6, 12, 16 though with a tiny minimum between 



the VBS and SF regions. It appears as if physical order 
parameters belong to some four dimensional surface and 
may diffuse over it without large effective action barriers. 

The minimum gets more pronounced for L = 24, 32, 
and breaks the distribution into well separated peaks 
for L > 48, i.e. the energy barrier flnally gets large 
enough to localize the order parameter in one of the 
phases. The two-peak structure is a smoking gun evi- 
dence that the transition is ultimately weak flrst-order. 
Further evidence is provided by finite-size scaling of the 
critical-point distribution functions, i.e the dependence 
of Ps{B — 0,L) (normalized to Ps{D = 0) as in Fig. ^ 
and the position of the P{B) distribution maximum, 
Bmax{L). The data in Fig.Osuggest that both quantities 
saturate to finite values in the thermodynamic limit. 

The study of the SF-CB transition along the K = 0.7, 
D = Q line reveals a remarkable similarity with the 
SF-VBS point. At J^ = 0.4184(2) the superfluid and 
solid orders exchange places with pronounced fluctua- 
tions of both order parameters in small systems. These 
fluctuations are almost identical to what is expected 
for the 0(3)-symmetric Heisenberg point where the dis- 
tribution function for the staggered order parameter, 
P{M < M,„ax) — P{Sz) is a step function. However, in 
larger systems a minimum in P{M) is developed. The 
double-peak structure of P{M) for L = 32 leaves no 
doubt that we are actually dealing with the weak first- 
order transition. Apparently, system sizes in the study 
of the SF-CB point in Ref. |j were too small to see the 
first-order transition. 

At certain conditions (which with J-current model can 
be easily achieved by adding appropriate terms) the su- 
persolid (SS) phase may intervene between SF and S 
phases. In this case, the vicinity to a (quasi-)0(3)/0(4) 
symmetric point with broken symmetry may render the 
SF-SS and SS-S transitions also first order, while nor- 
mally one might expect them to be of the second order. 

In conclusion, we note that weak first-order transitions 
discussed here can hardly be an artifact of the J-current 
model since they reveal themselves on large space-time 
scales. 
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